Prediction of narrow N* and A* resonances with hidden charm above 4 GeV 
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The interaction between various charmed mesons and charmed baryons are studied within the 
framework of the coupled channel unitary approach with the local hidden gauge formalism. Several 
meson-baryon dynamically generated narrow N* and A* resonances with hidden charm are predicted 
with mass above 4 GeV and width smaller than 100 MeV. The predicted new resonances definitely 
cannot be accommodated by quark models with three constituent quarks and can be looked for at 
the forthcoming PANDA/FAIR experiments. 
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Up to now, all established baryons can be ascribed 
into 3-quark (qqq) configurations [lj, although some of 
them were suggested to be meson-baryon dynamically 
generated states 04a] or states with large (qqqqq) compo- 
nents [9l4llT| . A difficulty to pin down the nature of these 
baryon resonances is that the predicted states from vari- 
ous models are around the same energy region and there 
are always some adjustable ingredients in each model to 
fit the experimental data. In this letter, we report a study 
of the interactions between various charmed mesons and 
charmed baryons within the framework of the coupled 
channel unitary approach with the local hidden gauge 
formalism. Several meson-baryon dynamically generated 
narrow N* and A* resonances with hidden charm are pre- 
dicted with mass above 4 GeV and width smaller than 
100 MeV. The predicted new resonances can be looked 
for at the forthcoming PANDA/FAIR experiments [l3J |. 
If confirmed, they definitely cannot be accommodated by 
quark models with three constituent quarks. 

We follow the recent approach of Ref. [12j and extend 
it from three flavors to four. We consider the PB — > 
PB and VB — > VB interaction by exchanging a vector 
meson, as shown by the Feynman diagrams in Fig. [T] 

The effective Lagrangians for the interactions involved 
are 0: 

Cvvv = ig{V»[V v ,d,y v ]) 
Cppv = -ig{V^[P,d^P\) 
Cbbv = g((B^[V^B]) + (B^B)(vn) (1) 

where P and V stand for pseudoscalar and vector mesons 
of the 16-plet of SU(4), respectively. 

Under the low energy approximation, the three mo- 
mentum versus the mass of the meson can be neglected. 
We can just take the 7° component of Eq. ([T]). The three- 
momentum and energy of the exchanged vector are both 
much smaller than its mass, so its propagator is approx- 
imately g^/My. Then with g = My/2/ the transition 
potential corresponding to the diagrams of Fig. [1] are 
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FIG. 1: The Feynman diagrams of pseudoscalar-baryon (a) 
or vector- baryon (b) interaction via the exchange of a vector 
meson. Pi, P 2 is D~ , D° or D~ , and Vi, V 2 is D*~ , D*° or 
D* s ~ , and Si, B2 is E c , Aj, S c , Hj, or Q, c , and V* is p, K* , <f> 

or lo. 



given by 
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ab(V 1 B 1 ->V 2 B 2 ) 
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= ^(En+Evjei-Z, (3) 



where the a, b stand for different channels of Pi(V\)Bi 
and Pity'ijBi, respectively. The E is the energy of cor- 
responding particle. The e is the polarization vector of 
the initial or final vector. And the e\ 2 component is ne- 
glected consistently with taking p/My ~ 0, with p the 
momentum of the vector meson. The C a b coefficients 
can be obtained by the SU(4) Clebsch Gordan Coeffi- 
cients which we take from Ref. [14|. We list the values 
of the C a b coefficients for PB — s- PB with isospin and 
strangeness (I, S) = (1/2, 0) and (0, -1) explicitly in Ta- 
ble I and Table II, respectively. 

With the transition potential, the coupled-channel 
scattering matrix can be obtained by solving the coupled- 
channel Bethe-Salpeterequation in the on-shell factoriza- 
tion approach of Refs 



r equ 



T= [l-VG]- l V 



(4) 



with G being the loop function of a meson (P), or a vec- 
tor (V), and a baryon (B). The poles in the T matrix are 



looked for in the complex plane of \fs. The e*i • e*2 factor 
of Eq. ([3]) factorizes out also in T. Those appearing in 
the first Riemann sheet below threshold are considered 
as bound states whereas those located in the second Rie- 
mann sheet and above the threshold of some channel are 
identified as resonances. 

TABLE I: Coefficients C ab in Eq. © for (I, S) = (1/2,0) 
DT, C DA+ g c N tyN V N -q'N KT, KA 
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For the G loop function, there are usually two ways to 
regularize it. In the dimensional regularization scheme 
one has [g, Il2j 
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where g is the four-momentum of the meson, P the total 
momentum of the meson and the baryon, s = P 2 and 
q denotes the three momentum of the meson or baryon 
in the center of mass frame, fi is a regularization scale, 
which we take to be 1000 MeV here. Changes in the 
scale are reabsorbed in the subtraction constant a M to 
make results scale independent. 

The second way for regularization is by putting a cutoff 
in the three-momentum. The formula is I3fl: 



the same close to threshold and they take the same value 
at threshold. This limits the a M to be around -2.3 with 
the corresponding A around 0.8 GeV, values which are 
within the natural range for effective theories [S(. Since 
varying the G function in a reasonable range does not 
influence our conclusion qualitatively, we present our nu- 
merical results in the dimensional regularization scheme 
with a M = —2.3 in this letter. 

From the T matrix for the PB -► PB and VB -> VB 
coupled-channel systems, we can find the pole positions 
zr. Six poles are found in the real axes below threshold 
and therefore they are bound states. For these cases the 
coupling constants are obtained from the amplitudes in 
the real axis. These amplitudes behave close to the pole 



T a b = 
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We can use the residue of T aa to determine the value of 
g a , except for a global phase. Then, the other couplings 
are derived from 



r ,9a^ab\ 

lim (— — ) 

\/B^rz R ±aa 



(8) 



The obtained pole positions zr and coupling constants 
g a are listed in Tables IIIII and IIVI Among six states, 
four of them couple only to one channel while two states 
couple to two channels. 
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4269 


2.85 


DA+ 





(0,-1) 


4213 
4403 


D 3 A+ 

1.37 




DE C 

3.25 




DE' C 


2.64 


TABLE III: Pole positions 
the states from PB — > PB. 


zr and coupling constants g a for 
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4418 


ZT£ C 
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D'At 





(0,-1) 


4370 
4550 


D* s At 

1.23 




D*E C 

3.14 




D*E' C 


2.53 



TABLE IV: Pole position and coupling constants for the 
bound states from VB — > VB. 
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where lu p = yjq 2 + M'p, lo b = \J a 2 + M/f, and A is the 
cutoff parameter in the three-momentum of the function 
loop. 

For these two types of G function, the free parameters 
are a M in Eq. ([5]) and A in Eq. ©. We choose a M or 
A so that the shapes of these two functions are almost 



As all the states that we find have zero width, we 
should take into account some decay mechanisms. Thus, 
we consider the decay of the states to a light baryon 
plus either a light meson or a charmonium through heavy 
charmed meson exc hang es by means of box diagrams as 
it was done in [15|, ll6(. Coupling to these additional 



channels with thresholds lower than the masses of pre- 
viously obtained bound states provides decay widths to 



these states and modifies the masses of these states only 
slightly. The results are given in Tables IVlandlVIl We do 
not consider the transitions between VB and PB states 
because in our t-channel vector meson exchange model 
they involve an anomalous VVP vertex which is found to 
be very small [15| • The transitions between VB and PB 
states through t-channel pseudoscalar meson exchanges 
are also found to be very small. As an example, we esti- 
mate the partial decay width of our DS C molecular state 
iV*g"(4265) to the J/tp p final state through the t-channel 
pseudoscalar D° meson exchange as shown by Fig. [3] 
Following a similar approach as in Ref. 17], the par- 



tial decay width is about 0.01 MeV, which is 3 orders of 
magnitude smaller than the corresponding decay to rj c p 
of 23.4 MeV. 



(i,s) 


m r 




Ti 


(1/2,0) 


nN 
4261 56.9 3.8 


7]N 

8.1 


n'N KT, rj c N 
3.9 17.0 23.4 


(o,-i) 


KN 
4209 32.4 15.8 
4394 43.3 


ttE 
2.9 
10.6 


rjA rj'A KE -q c A 
3.2 1.7 2.4 5.8 
7.1 3.3 5.8 16.3 


TABLE V 
width (I\) 


: Mass (M), total width (F), and the partial decay 
for the states from PB — >• PB, with units in MeV. 


(I, s) 


M r 




r, 


(1/2,0) 


pN 
4412 47.3 3.2 


10.4 


K*E J/tpN 
13.7 19.2 


(o,-i) 


K*N 
4368 28.0 13.9 
4544 36.6 


p£ 
3.1 

8.8 


wA cpA K*E J/xpA 
0.3 4.0 1.8 5.4 
9.1 5.0 13.8 



TABLE VI: Mass (M), total width (r), and the partial decay 
width (Ti) for the states from VB — >• VB with units in MeV. 
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FIG. 2: Feynman diagram for AT* 8 + (4265) -> J/tp p. 



In order to look for these predicted new V* and A* 
states, we estimate the production cross section of these 
resonances at FAIR. With a p beam of 15 GeV/c one has 
^/s = 5470 MeV, which allows one to observe N* reso- 
nances in pX production up to a mass Mx — 4538 MeV 
or Y* hyperon resonances in AY production up to a 
mass M Y ^ 4355 MeV. We take V c *+(4265) as an 
example. Its largest decay channel is r\ c p. Following 
the approach as in Ref. [181], we calculate its contribu- 
tion to pp — > pprjc through processes as shown in Fig. 
[3Ja,b) and also those from the conventional mechanism 
as shown in Fig. [3jc,d). For the conventional mecha- 
nism, the pprj c couplingis determined from the partial 
decay width of r\ c — > pp [1[ . For the new mechanism with 
the iV* e + (4265), its couplings to r\ c p and ixp are deter- 
mined from its corresponding partial decay widths listed 
in Table [Vj It is found that while the conventional mech- 
anism gives a cross section about O.lnb, the new mecha- 
nism with the V*^ (4265) results in a cross section about 
O.lfib, about 3 orders of magnitude larger. With the de- 
signed luminosity of about 10 cms for the p beam 
at FAIR [l3j, this corresponds to an event production 
rate of more than 80000 per day. With branching ratios 
for rj c — > KKir, r/mr, K + K~tt + tt~ , 2tt + 2tt^ of a few 
percent for each channel, the iV* E + (4265) should be easily 
observed from the rj c p and r\ c p invariant mass spectra for 
the pp — > pprjc reaction by the designed PANDA detec- 
tor [l3J]. The iV c * E + (4265) should also be easily observed 
in the pp — > ppJ ftp reaction with clean J/tp signal from 
its large decay ratio to e + e~ and /i" 1 " pT although the pro- 
duction rate is about 3 orders of magnitude smaller than 
the pp — > pprjc process. 
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It is very interesting that the six TV* and A* states are 
all above 4200 MeV, but with quite small decay widths. 
In principle, one might think that the width of these 
massive objects should be large because there are many 
channels open and there is much phase space for decay. 
However, because the hidden cc components involved in 
these states, all decays within our model are tied to the 
necessity of the exchange of a heavy charmed vector me- 
son and hence are suppressed. If these predicted nar- 
row N* and A* resonances with hidden charm are found, 
they definitely cannot be accommodated by quark mod- 
els with three constituent quarks. 



FIG. 3: Feynman diagrams of the reaction pp — > ppr\ c 



The D*T. C molecular state V*(4415) has a large decay 
branching ratio to J/V^ P- Its contribution to the pp — ► 
ppJ/ip reaction is estimated to be around 2nb, about one 
order of magnitude larger than the contribution from the 
V* e + (4265), and hence should be observed more clearly 
in this reaction. Similarly, the predicted DjAf-DE c 
coupled-channel bound state A* e (4210) states should 
be easily observed in pp — > AArjc reaction at FAIR. 
The other three predicted A* E resonances have too high 



masses to be produced at FAIR, but may be studied in 
some future facilities with higher p beam energies by the 
pp — > AAr/ c or pp — > AAj/ip reactions. This is an ad- 
vantage for their experimental searches, compared with 
those baryons with hidden charms below rj c N threshold 
proposed by other approaches [19J ] . 

In summary, we find two N* 5 states and four A* g 
states from PB and VB channels by using the local hid- 
den gauge Lagrangian in combination with unitary tech- 
niques in coupled channels. All of these states have large 
cc components, so their masses are all larger than 4200 
MeV. The width of these states decaying to light meson 
and baryon channels without cc components are all very 
small. On the other hand, the cc meson - light baryon 
channels are also considered to contribute to the width 
to these states. Then rj c N and r; c A are added to the 
PB channels, while J/ipN and J/i/jA are added in the 
VB channels. The widths to these channels are not neg- 
ligible, in spite of the small phase space for the decay, 
because the exchange D* or D* mesons were less off- 
shell than the corresponding one in the decay to light 
meson - light baryon channels. The total width of these 
states are still very small. We made some estimates of 
cross sections for production of these resonances at the 
upcoming FAIR facility. The cross section of the reaction 
pp — > pprjc and pp — > ppJ/ip are about 0.1//& and 0.2n6, 
in which the main contribution comes from the predicted 
A c * e (4265) and iV* g (4415) states, respectively. With this 
theoretical results, one can estimate over 80000 and 1700 
events per day at the PANDA/FAIR facility. Similar 
event rate is expected for the predicted A* e (4210) state 
in the pp — > AArj c reaction. These 3 predicted new nar- 
row N* and A* resonances should be easily observed by 
the PANDA/FAIR. The other 3 predicted A* g resonances 
will remain for other future facilities to discover. 

Although in the scheme of dynamical generated states 
these new N* e and A* E states are simply brothers or sis- 
ters of the well-known A* (1535) and A* (1405) in the 
hidden charm sector, their discovery will be extremely 
important. While for the A* (1535), A* (1405) and many 
other proposed dynamical generated states cannot clearly 
distinguish them from those generated states in various 
quenched quark models with qqq for baryon states and 
qq for meson states due to many tunable model ingre- 
dients, these new narrow N* and A* resonances with 
mass above 4.2 GeV definitely cannot be accommodated 
by the conventional 3q quark models, although how to 
distinguish these meson-baryon dynamically generated 
states from possible five-quark states needs more detailed 
scrutiny. The existence of these new resonances with hid- 
den charm may also have important implications to the 
long-standing puzzles relevant to charmonium produc- 
tion in various collisions involving nucleon in the initial 
state, such as the strikingly large spin-spin correlation 
observed i n pp elastic scattering near charm production 
threshold [20] and difficulties in reproducing the cross 



sections and polarization observables of JM production 
from high energy pp, pp and jp reactions |2ll.l22|. These 
issues deserve further exploration. 
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